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^^ ■ Gravitons produced from quantum vacuum fluctuations during an inflationary stage in the early 

f^ ' Universe have zero entropy as far as they reflect the time evolution (squeezing) of a pure state, 

C^ I their large occupation number notwithstanding. A non-zero entropy of the gravitons (classical 

^SJ . gravitational waves (GW) after decoherence) can be obtained through coarse graining. The latter has 

, , ' to be physically justified and should not contradict observational constraints. We propose two ways 

O , ' of coarse graining for which the fixed temporal phase of each Fourier mode of the GW background 

still remains observable: one based on quantum entanglement, and another one following from the 

presence of a secondary GW background. The proposals are shown to be mutually consistent. 

OO ' They lead to the result that the entropy of the primordial GW background is significantly smaller 

, than it was thought earlier. The difference can be ascribed to the information about the regular 

O^ ' (inflationary) initial state of the Universe which is stored in this background and which reveals 

^ ' itself, in particular, in the appearance of primordial peaks (acoustic peaks in the case of scalar 

V^ perturbations) in the multipole spectra of the CMB temperature anisotropy and polarization. 

^ . 

5 ! I. INTRODUCTION 

On ' What is the entropy of primordial fluctuations from which all inhomogeneities in our Universe are assumed to 

f^ ' originate? This fundamental problem becomes especially intriguing in the framework of the inflationary paradigm. 
^^. Indeed, in that case primordial fluctuations are generated from vacuum quantum fluctuations. In this work, we shall 
J ' concentrate on a primordial gravitational-wave (GW) background in inflationary cosmological models, although a 
(2J[), similar discussion can be applied to scalar perturbations in such models as well. It is well known that in this case the 
|V ' unitary evolution of the quantum state of the fluctuations outside the Hubble radius leads to an effective quantum-to- 
. ^H ', classical transition of a very specific kind: the fluctuations become indistinguishable from classical fluctuations with 
S^ stochastic amplitude and fixed temporal phase once the so-called decaying mode of the perturbations is neglected [Q . 

5-H The approximation of simply omitting the decaying mode and then introducing the equivalent classical stochastic 

fluctuation is sufficient for calculating correctly the amplitudes of observable inhomogeneities in the Universe. It 
does not require any explicit account of interaction with an environment. In particular, inflationary predictions for 
the amplitude, statistics and quasi-classical temporal behaviour of perturbations are independent of this interaction 
(at least, for sufficiently large scales). Note that the possibility to get a decoherence-independcnt prediction for 
observable quantities (inhomogeneities) is a specific property of the very mechanism of generation of inhomogeneities 
in the inflationary scenario of the early Universe, and does not occur for general quantum systems, as has been shown 
in ||] with the help of the Wigner function and in [g[ Q by a discussion of thought experiments involving slits. Also, 
it may be attributed to the fact H] that amplitudes of perturbations constitute an (almost ideal) pointer basis in 
this case. Finally, this approximation fits well into the consistent-histories framework: as the decaying mode becomes 
negligible, probabilities can be consistently assigned to classical trajectories ("histories") of perturbation modes in 
phase space |^ . 

However, this does not mean that this approximation is sufficient for the calculation of all quantities. In the present 
paper we consider one important (though not directly measurable) quantity for which it is definitely insufficient and 
for which one has to consider actual mechanisms of environment-induced decoherence - the entropy of cosmological 
inhomogeneities. The entropies of an initial quantum perturbation in a pure strongly squeezed state and of its 
(approximate) equivalent (a classical perturbation with a stochastic amplitude, but a fixed phase) are in fact both zero. 
It is, however, clear that, as a result of an effective coarse graining produced by environment-induced decoherence, a 
pure quantum state becomes mixed, while the fixed classical phase of the classical system acquires some small stochastic 
part, so a non-zero entropy should arise in both quantum and effective classical descriptions of perturbations. 
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The main question, then, is: what is the correct coarse graining? Any coarse graining introduced by hand must be 
justified by some concrete physical mechanism. Also, it must be in agreement with present and future observational 
constraints. Indeed, as was already emphasized in Q, the coarse graining that one chooses for scalar (adiabatic) 
perturbations can be crucial with respect to the appearance, or not, of periodic oscillations - acoustic peaks - in the 
CMB anisotropy multipoles Cf . The same refers to the primordial GW background where the corresponding periodic 
oscillations (with a different period) are more easily seen in the CMB polarization multipoles Cpi [|l| Q (in this case, 
it is more appropriate to call them primordial peaks). 

In earlier papers [g] where the entropy of cosmological perturbations was considered, an averaging procedure was 
performed with respect to the squeezing angle ipk of each Fourier mode, leading to the entropy Sk = 2r/c per mode, 
where r^, is the squeezing parameter. Such a "high" entropy would, however, spoil the generic prediction of the 
inflationary scenario about the dominance of the quasi-isotropic (or growing, in the case of scalar perturbations) 
mode at the moment of the last Hubble-radius crossing, which represents the equivalent classical expression of the 
remaining coherence between the field-mode amplitude and its momentum. As was shown in pi , the latter coherence 
persists long after the loss of quantum correlation between quasi-isotropic and decaying modes of perturbations (and 
may survive even up to the present moment for sufhciently large perturbation wavelengths). So, in the case of 
cosmological perturbations generated during inflation, the quantum decoherence sufficient for the classical description 
of the perturbations occurs long before the time of complete relaxation (if this relaxation occurs at all) . 

It is just this property that leads to the appearance of the above-mentioned peaks in the CMB fluctuation spectra. 
The advent of high-precision cosmological observations in general, and accurate measurements of the CMB fluctuations 
in particular, has dramatic consequences on this problem, which seemed remote only a few years ago: the existence 
of peaks is incompatible with complete randomization of the temporal phase of the cosmological fluctuations and 
severely restricts the choice of possible coarse grainings. 

Since, as is well known, the equations for small scalar (adiabatic) perturbations and GW (tensor) perturbations 
superimposed on a Friedmann-Robertson- Walker (FRW) background can be reduced to the equation of a real mini- 
mally coupled scalar field 0(r, t), massless in the case of GW, while there is a mass term and a sound velocity in the 
case of adiabatic perturbations, we shall further consider this auxiliary field. The FRW background is assumed to be 
spatially fiat: 

ds^ =dt^ ~a^{t)df , (1) 

where dt^ is the three-dimensional Euclidean interval (c = ?i = 1 is assumed throughout the paper). The quantum 
field (j) can then be decomposed into Fourier harmonics 

0(k) = (0(k,i)a(k)+r(-k,i)at(-k))e'''---0t(_k) ^ (2) 

with a(k) and a^(k) being the usual annihilation and creation operators in the Fock space, respectively. 

Let us take a two-mode system (/)(k), 0(— k). This system exhibits a semiclassical behaviour by itself soon after 
the first Hubble-radius crossing during an inflationary stage due to the disappearance of the decaying mode (in the 
Heisenberg representation), or equivalently due to large squeezing (in the Schrodinger representation). More precisely, 
once the wavelength of the mode A = a{t)/k, k = |k| is much bigger than the Hubble radius Rh = H{t)~^, H = a/ a, 
we can make (/)(k, t) = (f>{k,t) real by a time-independent phase rotation. Then, after the omission of the decaying 
mode, the quantum mode operator (g) becomes equivalent to the Gaussian stochastic quantity [y 

0(k) = <^*(-k)=c(k)</>(fc,t)e^•'^ (3) 

where c(k) — c*(— k) is a complex Gaussian stochastic fleld with the following non-zero correlations: 

(c(k)c*(k')) = <5(k - k') , (c(k)c(k')) = S{k + k') . (4) 

The stochastic field (^ has still zero entropy since its volume in phase space is equal to zero. However, interactions 
with other fields, even tiny ones, are unavoidable. This is a first physical process which can lead to coarse graining. It 
occurs on very short decoherence timescales, rendering the field-amplitude basis a classical "pointer basis" |§, Q - the 
robust basis with respect to the environment j^, |l^. Once A ^ Rh, the dynamical infiuence of such environmental 
fields should be negligible since causal processes are prohibited. In particular, the quasi-isotropic mode is not affected. 
However, for decoherence it is just sufficient to destroy the quantum correlation between the decaying and quasi- 
isotropic parts of perturbations. Since the decaying mode is "local" (in the sense that it may be affected by local 
processes), decoherence may be easily achieved by coarse graining of quantum entanglement between (the decaying 
mode of) perturbations and the environment. As was pointed out above, the relaxation timescale exceeds by far 
the decoherence timescale (the latter being of the order of H~^ at the first Hubble radius crossing M). Since the 



maximuiii entropy 2rfe is associated with the relaxation timescale, the actual entropy arising from entanglement, 
although non- vanishing, should be much smaller. 

In the following we shall consider two different coarse grainings which are physically relevant. First, a realistic 
situation is considered where interaction with environmental fields is efficient in the suppression of off-diagonal terms 
in the density matrix, but docs not produce any significant back reaction (Sec. II). Since it appears that details 
of the interaction arc not relevant, we model its effect through the introduction of a phenomenological parameter 
suppressing the off-diagonal terms. Second, the entropy growth due to the loss of information about the primordial 
GW background which occurs as a result of the appearance of a secondary GW background after the second Hubble 
radius crossing is calculated in Sec. III. This effect arises because we cannot distinguish one kind of gravitons from 
another at the present epoch. So, in the first case a non-zero entropy arises because the primordial GW background is 
directly affected by an environment, while in the second case it is simply "polluted" by a secondary GW background 
emitted by some environment. It will be shown that both coarse grainings yield in general Sk significantly smaller 
than 2rk- Sec. IV contains conclusions and discussion, including the relation of this problem to the problems of 
entropy growth in open systems and in the Universe. 

II. ENTROPY DUE TO COARSE GRAINING OF QUANTUM ENTANGLEMENT 

Let us consider a squeezed vacuum state which describes the behaviour of primordial fiuctuations in the absence of 
interactions with the environment. It has the form 

V'feo=(^j exp{-[nn + ini]\yk\^) , (5) 

where yk is a Fourier component of the rescaled quantum field y(r, t) — a(i)0(r, t). Since modes with different wave 
vectors k decouple, we shall sometimes skip the index k (or k) in the following. The wave function (ra) can be written 
in terms of the squeezing parameters r and ip in the form 

_/ 2k y/^ / 1-e^^^tanhr 2\ 

^° ~ V7r[cosh2r-Hcos2.^sinh2r]y exp I - j^ ^_ g2i./> ^anh r i ' 

Squeezing can equivalently be expressed in terms of "particle creation" with average particle number iV(k) « e^''/4. 
We are interested in the regime after the first Hubble-radius crossing (occurring during an inflationary stage) when 
A^(k) > 1. 

To calculate the entropy it is sufficient to consider one of the two modes contained in the complex amplitude 
?/(k, i) = ?/*(— k, i) (in the following, y is a real function). The entropy is additive for the various modes. The density 
matrix corresponding to the pure state (ra) is given in the field-amplitude representation by 



PQ[y,y) = y-^cxp i — — 



-'-^iy-y'f--^^i{y-y'){y + y')-^{y + y'f] ■ (7) 



The coefficient in front of {y — y'Y is a measure of quantum coherence (size of non-diagonal elements), while the 
coefficient in front of {y + y'Y is related to (Ay)^ and measures the extension in configuration space of a fictitious 
ensemble described by the density matrix (details can be found in Appendix A2.3 of pO|). Since for large squeezing 
U,ji becomes very small, both coherence and extension in phase space become large. This coherence cannot be 
distinguished from a classical random process, see for instance the thought experiment involving slits in P, H, p|. 
Coupling to environmental degrees of freedom |^, Q will decrease the coherence length in j/, i.e. increase the coefficient 
in front of {y — y'Y ■ Although details may be very complicated, a very good approximation is to multiply po by a 
Gaussian factor which suppresses off-diagonal elements but leaves diagonal elements (probabilities) untouched. This 
corresponds to the situation where dynamical back reaction is small ( "ideal measurement" ) . Then, instead of (0) , the 
following density matrix is obtained: 



pdy^y') ^ po [y, y') exp f - - (2/ - y')^ j 



(8) 



where ^ ^ f7/j. For large squeezing, this condition reads 



^ »1. (9) 



In the following, it is referred to as the decoherence condition. 

To study correlations that might be present between y and (the Fourier transform of) the momentum p, it is 
preferable to calculate the Wigner function. This is found to be 

VK,(„)4»p(-2<|0|!-2n«.^)^ (10) 

The variances Ay and A(p — p^;) can be found from the extension of the corresponding contour ellipse which describes 
the extension of the (apparent) ensemble. This ellipse becomes a circle in the vacuum case (r = 0, ^ = 0), if we 
use the axes p/k, y. It will be convenient to introduce the variables y = 2\/ky, p = -yrP, and we denote the major 
half axis of the Wigner ellipse in the y, p plane by a and the minor half axis by /?. In the absence of environment 
(f — 0), one has ao = e^, Po — e~'^ O. For ^ 7^ one generally finds rather complicated expressions, from which 
simple expressions are recovered in the limit of large squeezing. One finds for e^'' — > c» 

a « e^ /? « yi » /3o . (11) 

While the size of the major half axis remains the same, the size of the minor half axis becomes much bigger than 
before due to decoherence. 

To preserve a correlation between y and p one has to demand that a remains much bigger than /3, or: 

^ < 1 . (12) 



ke^ 



We shall refer to (12) as the correlation condition. The surface of the ellipse (divided by n) is A — a/3 = e''-\/^/A:. 



Thus, one would expect the following relation between entropy and volume in phase space: 

S^lnA^ (1/2) ln(e2'-,f /fc) . (13) 

This is what we shall demonstrate in the following. 

The von Neumann entropy connected with the density matrix (ph is given by 

S=-Tiip^\np^). (14) 

It was calculated in |1^ for an arbitrary Gaussian density matrix, see also Appendix A2.3 in ||l^. The result is 

S — —\npo Inq , (15) 



Po 



where 



2V^. V^R + C-V^ ,^„. 

Pn = — , : , q = — , ; . lo 

Using e^ ^ CO and ^ ^ f2_R, one finds 

1 e'^^'E 1 N£ 

5«l-ln2 + -ln^ = l + -ln^, (17) 

in accordance with the expectation S* ~ In A. Applying the decoherence condition (0), one finds 

5> l-ln2«0.31 , (18) 

where ^ holds here in a logarithmic sense (it directly holds for the number of states e"^). Note that this lower bound 
on the entropy corresponds to the loss of less than one bit of information. This is consistent with previously known 
results on decoherence in quantum mechanics. For example, in recent quantum-optical experiments p2[ , decoherence 
starts if, on average, one photon is lost. Thus, it may be expected that a minimal entropy Smin ~ hi 2 per mode 
would be sufficient to guarantee decoherence in the present case, too. 
Applying the correlation condition (|l2|), we get 

S < 2r, (19) 



where again ^ holds in a logarithmic sense. It is evident that the entropy must be much smaller than the maximum 
value 2r which is found by integrating over the squeezing angle |g[| (and which would just mean a = /3 « e'' for the 
Wigner ellipse). 

Note that, though our initial formula S = \nA (proposed earlier in ||6[, too) is the same as the one introduced 
by Rothman and Anninos p3[ , our final result (Q) is drastically different from the conclusion of the recent paper 
by Rothman ||lj| that this formula for the entropy leads to results identical to those obtained in [g|. The reason for 
this difference is evidently the fact that in [O, nj] , the total phase-space volume of a state with a given energy was 
calculated. This does not properly account for the squeezed nature of the Wigner ellipse in our case. 

Although details about ^ may be complicated, it is natural to expect that the coherence length ^^^/^ is not smaller 
than the width of the ground state for r = 0, at least during the inflationary stage, so that the quantum state really 
remains squeezed in some direction as compared to the ground vacuum state. This means that ^ < k and consequently 
S < r. Therefore, pure decoherence (without dynamical influence) can never totally smear out the Wigner ellipse to 
achieve S ~ Smax = 2r. For complete randomization, one would thus have to invoke, for example, a thermal bath at 
a sufficiently high temperature. Such a model was discussed in ||l^], and their results seem to be consistent with our 
general treatment. Conditions similar to (y) and (|l2| ) have been frequently discussed in quantum mechanics, see e.g. 
Eqs. (6.32)-(6.36) in |l|. 

III. ENTROPY DUE TO A SECONDARY GW BACKGROUND 

Let us consider now a different coarse graining due to the presence of a secondary GW background which is 
generated by different matter sources in a causal way after a given mode of the fluctuations has re-entered the Hubble 
radius |0, Qq]. For primordial GW, as shown in |H], the coherence with respect to the squeezing angle is maintained for 
a considerable time even after the second Hubble-radius crossing. On the other hand, the secondary GW background 
is expected to have a uniformly distributed phase. It is natural, and rather general, to assume that a k-mode of 
this secondary background is described by a density matrix p'^(k) which is diagonal in the occupation number basis, 
p*' = ^^ Wn\n){n\ (here the argument k is omitted, and n = 0, 1...). This means in particular complete randomization 
of the temporal phase, or the absence of a preferred direction in phase space. As gravitons of the primary (squeezed) 
background are indistinguishable from those belonging to the secondary background, we expect some information loss 
when both backgrounds are mixed. The mean occupation number 

oo 

n(k) = ^ nw;„ , (20) 

though expected to be significantly lower than that of theprimary background iV(k), may nevertheless be large, too. 
We shall show that, loosely speaking, p* corresponds to (g) with 

C :^ i^fcn(k) . (21) 

Let us estimate the corresponding entropy. We consider the typical volume F in (half of) phase space occupied 
by the system. For the isolated primary background, this volume Fo is the minimal one which corresponds to zero 
entropy. The secondary background occupies a much larger volume which corresponds to circles in phase space with 
the approximate radius \/{y^ y'^^), where 

k {y^ y^t) = „(k) + 1 . (22) 

Here the averaging process involves also time averaging in addition to quantum average. Therefore at any time, 
the width of F in the squeezed direction of the primary background will be given by the radius of the secondary 
background, while the elongated direction of the primary background remains dominated by the primary background. 
However, since the vacuum part 1/2 should not be counted twice (it has been already counted and "squeezed" in the 
primordial part of the mode), only newly created gravitons n(k) produce an additional noise in the initially squeezed 
direction in phase space. If A^(k)7i(k) ^ 1, then F, the volume in phase space occupied by the whole system, can be 
estimated as follows: 

T'^^al n{k) ^ iV(k)n(k) » T^ = ^ , (23) 

where Fq is the phase-space volume of the initial vacuum state. The corresponding entropy S'(k) can be estimated as 

5(k) = In — = rfc + ln2 + i Inn(k) ^ rk +ln2n - -\n— + -Infl'' , (24) 

Fo 2 2 Ec 2 



where we have used the physical quantity n'^luj) = — ^^ = uj^n^/n'^ec {uj — k/ao) which describes the secondary 
stochastic background, while e^ denotes the critical density. Comparing ( p4) ) with (p^), the factor 16/e^ is found in 
(pl|). It follows that the condition iV(k)n(k) ^ 1 used above corresponds to the decohcrence condition (pi), while the 
condition that the secondary background is much smaller than the primordial one, -/V(k) 3> '^(k), corresponds to the 



correlation condition (12) 



We see from (|24|) that S'(k) < 2rk as long as n(k) < iV(k), and also S'(k) > 5* w Inn(k) (for n(k) > 1). On 
large cosmological scales, we have N^k) ~ lO^"*', r^ ~ 115, and we have approximately for most models -/V(k) ex fc""* 
for 10^^^ <^ V — ui/2tt <^ 10^'^Hz. We expect that very little entropy is produced by this mechanism on large scales, 
hence S'(k) ^ 2rk for these scales. 
Several possibilities can be distinguished now. 

1) n(k) < 1. 

This case, certainly plausible on very large cosmological scales, yields S'(k) <C r^. 

2) n(k) « 1. 
Now we have 

S'(k) « rfc = -Sjnax • (25) 

Note, however, that n(k) <C A^(k) and we still have a significant squeezing in phase space. Therefore, the standing- 
wave behaviour of the primordial GW still has definite observational consequences. 

3) 1 < n(k) < A^(k). 

In that case we have r^ < S'(k) < 2rfc. However, like for the preceding case, the system remains highly squeezed. 
Therefore, with respect to observations, this case is very similar to the preceding ones. 

4) n(k) > N{k). 

Now, the squeezing of the primordial background is not apparent anymore in observations. When n(k) ^ N(\i.), the 
squeezing is completely "washed out" . 

So, we have here another concrete coarse graining of the primordial GW stochastic background. For all cases for 
which 1 <ti n(k) <C -^(k), which is a reasonable assumption on large cosmological scales, the resulting entropy is thus 
significantly smaller than 2rfc. The crucial point is that a value r^ < S(k) < 2rk is still fundamentally different from 
S'(k) = Smax = 2rfc. In the first case, the composite system of both backgrounds can still reflect the squeezing of the 
primordial GW background, while it does not in the second case. The peaks in the B-mode polarization of the CMB 
produced by the total GW background would be absent in the second case. 

IV. CONCLUSIONS AND DISCUSSION 

We have shown that both methods of coarse graining proposed yield for the entropy per mode the result S'(k) <C 2rk 
(in the logarithmic sense), which differs strongly from the results obtained earlier M. The reason is that our coarse 
grainings do not destroy the classical correlation (the almost deterministic temporal phase of the GW and of the 
adiabatic perturbations) remaining after quantum decohcrence, which, as physical estimates show [3[, persists for a 
rather long time (very long in the case of GW) after the perturbations have re-entered the Hubble radius. Thus, 
though the maximal entropy per mode S'(k) = 2rk might eventually be reached after complete relaxation, it is not 
reached at the recombination time for sufhciently long-wave scalar perturbations and even at the present time for 
sufficiently long-wave GW. As was mentioned above, this leads to observable effects: periodic peaks in the multipole 
power spectra of the AT/T-anisotropy and polarization of the CMB. Since the difference between entropies arising 
after two coarse grainings of the same system may be interpreted as the difference of amount of information loss due 
to these coarse grainings, the difference between previous results for the entropy of primordial GW background and 
our result may be interpreted as the amount of information contained in primordial peaks (acoustic peaks in the case 
of scalar perturbations). 

Note that the energy density of the primordial GW background calculated either with the quantum operators (g) 
or by the use of equivalent classical stochastic fields (S) is one half of the energy density of the GW background 
described by a density matrix which is isotropic in phase space, with the same average number of gravitons given by 
ri(k) (this, of course, was taken into account in all correct calculations of the energy density of GW produced during 
inflation, beginning with |lj]). Now for the entropy, S = Smax/'^ is also a distinguished, "average" case. However, we 
see that large deviations from Smax/'^ on both sides are possible; the entropy per mode can be both much larger or 
much smaller than r^ (in the logarithmic sense). 

Let us finally discuss the temporal growth of the entropy of cosmological perturbations, which is at least partly 
responsible for the arrow of time in the Universe [^. After the perturbations re-enter the Hubble radius, A^(k) 
remains constant (neglecting very small graviton-graviton scattering events and graviton absorption by matter), so 



the growth of entropy is mainly due to the growth of n(k, t) due to irreversible emission of gravitons by matter. 
Due to the expansion of the Universe, n(k, i) typically grows only as a power of time, so (S'(k, i) oc t~^. The same 
refers to the earlier period for A ^ Rh, but after the end of inflation. (In inflation, one has for long wavelengths 
S ~ H K:i constant and therefore S ex t.) Then, there is no secondary GW background and we may use ( p7| ) only. 
Under any reasonable assumptions about the time dependence of ^, ^(k, t) grows logarithmically, so ^(k, t) oc t^^, too. 
Comparing this behaviour with those considered in the interesting discussion pl| of the paper [ p2[ , we see that the 
case of a two-mode subsystem of cosmological perturbations is in some sense intermediate between the two different 
cases discussed in [^ (decoherence due to a chaotic behaviour of the subsystem itself or decoherence caused by the 
environment). Though our subsystem is not chaotic but only classically unstable, the law for S{k,t) is the same as 
for chaotic dynamical systems. 
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